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Abstract 

We study the limits of the empirical distributions of the eigenvalues of two n by n 
matrices as n goes to infinity. The first one is the product of m i.i.d. (complex) Ginibre 
ensembles, and the second one is that of truncations of m independent Haar unitary 
matrices with sizes rij x rij for 1 < j < m. Assuming m depends on n, by using the 
special structures of the eigenvalues of the two matrices we developed, explicit limits 
of spectral distributions are obtained regardless of the speed of m compared to n. In 
particular, we show a rich feature of the limits for the second matrix as rij In’s vary. 

Some general results on arbitrary rotation-invariant determinantal point processes are 
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1 Introduction 


In this paper we will study the limiting spectral laws of two types of random matrices. 
They are in the form of Xi • • • X^, which is called a product ensemble. The first type is the 
product of n X n Ginibre ensembles, that is, Xj, 1 < j < m, are independent and identically 
distributed (i.i.d.) Ginibre ensembles; review that the n x n matrix Xi is referred to as 
a Ginibre ensemble if its v? entries are i.i.d. standard complex normal random variables. 
The second kind corresponds to that Xj, 1 < j < m, are independent n x n matrices, 
each of which is a truncation of an Haar-invariant unitary matrix. We do not assume these 
matrices are of the same size. 

We first give the structures of the eigenvalues of the two matrices (Lemmas 11.21 and 11.31) 
by using a theory of the determinantal point processes. It is found that the absolute values 
of the eigenvalues are the product of i.i.d. Gamma-distributed random variables and the 
product of i.i.d. Beta-distributed random variables, respectively. 

Using the theory, assuming m depends on n, we obtain the limiting distributions of the 
eigenvalues of Xi • • • X^ for both cases as n —)• oo regardless of the speed of m = rUn- As 
m does not depend on n for the first case or m = 1 for the second case, some knowledge 
about their limiting distributions are known. Here our results hold for any choice of m„. 
For the product of truncations of Haar unitary matrices with different sizes, the limiting 
distributions are very rich. 

The essential role in the derivation of our results is the determinantal point process 
{Zi, - ■ ■ ,Zn}- For the two product ensembles above, their kernels associated with the point 
process are rotation-invariant. We then study it and obtain a general theory in Section 11.31 
They may be useful in other occasions. 

Before stating the main results, we need the following notation. 

• Any function g{z) of complex variable z = x + iy should be interpreted as a bivariate 
function of {x,y): g{z) = g{x,y). 

• We write g(z) dz = g(x, y) dxdy for any measurable set A C C. 

• Unif{A) stands for the uniform distribution on a set A. 

• For a sequence of random probability measures {v,Vn', n > 1}, we write 

Vn V if P{vn converges weakly to u as n —)• oo)=l. (1-1) 

When u is a non-random probability measure generated by random variable X, we simply 
write Vn ^ X. For complex variables {Zj] 1 < j < n} mentioned above, we write 

Qj = arg(Zj) € [0,27r) such that Zj = \Zj\ ■ (1-2) 

for each j. Let Li, • • • , Tn be some given random variables, each of which may also rely on 
n. We omit the index n for each Yj for clarity. Given a sequence of measurable functions 
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(1.3) 


hn{r), n > 1, defined on [0,oo), set 

^ n 1 ^ 

i=i i=i 

which are the empirical or counting measures of two sets of random variables. In particular, 
if hn is linear, that is hn{r) = r/an, where {a^ > 0; n > 1} is a sequence of numbers, we 
give special notation of the empirical measure of Zj's accordingly for this case by 

^ n ^ n 

^n = -Y.^Z,la^ and K = -^SY,la^. (1.4) 

j=i i=i 

Review the notation in dni). The symbol ^ii (8) /i 2 represents the product measure 
of two measures fii and fi 2 - Our general result is given as follows. 

Theorem l Let ^{x) >0 be a measurable function defined on [0, oo). Assume the den¬ 
sity of {Zi, ■■ ■ ,Zn)eC'^ is proportional to Y[i<j^k<n\^j-^k\^'Y[j=i^i\^j\)- LetYi,--- ,Yn 

be independent r.v.’s such that the density of Yj is proportional to y‘^^~^ip{y)I{y > 0 ) for 
every 1 < j < n. If {hn} are measurable functions such that Vn n for some probability 
measure v, then pin h with pi = Unif[0,27r] (g) n . Taking hn{r) = r/on, the conclusion 
still holds if “{ptniVn, pi,n)” is replaced by “{pLn,r'n, pi*,n*)” where pi* is the distribution of 
Re*® with (0,R) having the law of Uni f[0,27r] ®v*. 

Next we apply this theorem to the two product ensembles. The tasks are to obtain 
an explicit expression of hn, the probability measure n and to check the weak convergence 

“l, V 1,” 

h'n ^ ^ • 

1.1 Product of Ginibre Ensembles 

Given integer m > 1. Assume Xi,--- ,Xm are i.i.d. nx n random matrices and the 
entries of Xi are i.i.d. with the standard complex normal distribution CX(0,1). Let 
Zi, - ■ ■ , Zn be the eigenvalues of the product YVjLi It is known that their joint density 
function is 

n 

C \Zj - Zl\^Y\_'^rn{\Zj\) (1.5) 

l<j</<n j=l 

where C is a normalizing constant and Wm{z) has a recursive formula given by wi{z) = 
exp(—|zp) and 

Wjn{z) =27r f Wm-i(-) exp(-r^) — 

Jo ^rJ r 

for all integer m > 2; see, e.g., Akemann and Burda (2012). The function Wm{z) also has 
a representation in terms of the so-called Meijer G-function; see the previous reference. 
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Through investigating the limit of the kernel of a determinantal point process, Burda et 
al. (2010) and Burda (2013) showed that the empirical distribution of 1 < j < 

in the sense of mean value, converges to a distribution with density —for \z\ < 1. 
Under the condition that the entries of Xi are i.i.d. random variables with a certain 
moment condition, Gotze and Tikhomirov (2010) prove the above result in the sense of 
mean value. Bordenave (2011), O’Rourke and Soshnikov (2011) and O’Rourke et al. (2014) 
further generalize this result to the almost sure convergence. Our result next gives a weak 
convergence of eigenvalues Zj’s by allowing m = nin, that is, m depends on re, and the 
result holds regardless of the speed of rre^ relative to re. Review m- 

Theorem 2 Let {rUn > l; U > 1} he an arbitrary sequence of integers. Define 

1 ” 

Then fin ^ Urei/([0, 27r) (g) [0,1]) as re —)• oo. 

Theorem [2] implies that the angle and the length of a randomly picked pair (0j, ^ iZj-p/™") 
are asymptotically independent. Take rUn = rn. By the continuous mapping theorem, the 
above conclusion implies that, with probability one, the empirical distribution of —^Zj = 
\Zj\ e*®.’, 1 < J < re, converges weakly to the distribution of where {Q,R) 

follows the law C/rei/([0, 27 r) x [0,1]). Easily, Z := € C has density 

for \z\ < 1. This yields the conclusion mentioned before Theorem [2j In particular, taking 
rUn = 1, we have 

1 " 

-'^^Zj/^'-^Unif{\z\ <1}. ( 1 . 6 ) 

This gives the classical circular law. In general, let a square matrix be filled with i.i.d. 
entries (not necessarily Gaussian random variables) of mean zero and variance one, then 
the empirical distribution of its eigenvalues Zfis is asymptotically the circular law in the 
sense of (II.6p . There are a lot of results in this aspect, see, for example, Girko (1984), Bai 
(1997) and Tao and Vu (2010). A detailed exposition is given by Bordenave and Ghafai 
( 2012 ). 

The proof of Theorem [2] is based on Theorem [1] and a recent result by Jiang and Qi 
(2015) in which the exact distributions of \Zj\ are shown to be the products of independent 
Gamma-distributed random variables (see Lemma ll.2p . 

1.2 Products of Truncated Unitary Matrices 

Let rre, re and {Ij] 1 < j < m} be positive integers. Set nj = lj + n for 1 < j < m. Suppose 
{Ujj 1 ^ j < w-} are independent Haar-invariant unitary matrices where Uj is nj x nj for 
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each j. Let Xj be the upper-left n x n sub-matrix of Uj. Consider the product matrix 
= X^Xm-i • • • Xi. We write this way instead of the product of the matrices in the 
reverse order is simply for brevity of notation below. 


The joint density function for the eigenvalues Zi, - ■ ■ , of X*-™^ is derived by Akemann 
et al. (2014): 

n 



I ^n) 

= C 

n 



( 1 . 7 ) 




l<j<fc<n 

i=i 



for all Zj's with maxi<j<m 

\^j\ 

< h 

where C = 

1 rjfn 

n\ 1 lj =0 11/=0 

\ and Wm'"' 

can 

be recursively obtained by 








■ 27r 

/ w 


■ Im-l) ,s.du 

s G [ 0 , 1 ), 

(1.8) 



Jo 

u ’ 


with initial w^\s) = {I/tt 

)(1 

-s^) 

^-1/(0 < s 

< 1). The function Wm’ 

can be 


expressed in terms of Meijer G-functions. One can see Appendix C from Akemann et al. 
(2014) for details. The density in (jl.7p for the case m = 1 is obtained by Zyczkowski and 
Sommers (2000). 


We will consider the limit of the empirical distribution oi Zi, - ■ ■ ,Zn. Assume m = nin 
and 1 < re < mini<j<mrej. For convenience, we assume re-i, • • • ,nm are functions of re and 
all limits are taken as re —)• oo unless otherwise specified. The liming spectral distribution 
actually depends on the limit of functions F)i(x)’s defined below. 

Let { 7 ^; re > 1} be a sequence of positive numbers. Define 


m 



Note that Fn{x) is continuous and strictly increasing on [0,1], T„(0) = 0 and Fn{l) = 1. 
We will assume that Fn{x) has a limit F{x) defined on (0,1) such that 

F{x) is continuous and strictly increasing over (0,1), 

lima;^o L’(a;) = 0 and lima;|i F(a;) = 1 (1-10) 

(we denote f*{x) := ■^F~^{x) for 0 < a; < 1 if the derivative exists), or a limit F{x) defined 
on (0,1] satisfying 

F{x) = 1 for X G (0,1]. (1-11) 

Recall the notation in (jl.ip and “[/rez/(A)” standing for the uniform distribution 
on a set A. Write z = re*®. Evidently, z G C and {9,r) G [0, 27 r) x [0, 00) are one-to-one 
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correspondent. Although the empirical distributions considered in the following are targeted 
as functions of ( 0 , r), we will characterize their limits in terms of complex distributions since 
the arc law and the circular law etc are easily understood. 


Theorem 3 Assume there exists a sequence of numbers {7n} with 7n > 1 such that 
lim„^oo Fn{x) = F{x), X G (0,1), for some function F{x) defined on (0,1). Let Zi, - ■ ■ , 
be the eigenvalues ofX.^^\ be as in (11.31) with hn{r) = {r"^and bn = OS 

(a) . If Ill.lO]) holds and f*{x) = ■^F~^{x) exists for x G (0,1), then fin {Q,R) such 

that Z := iZe*® has density 0 < | 2 :| < 1. 

(b) . If m.ll]) holds, then fin {Q,R) such that Z = i?e®® has the law Unif{\z\ = 1}. 

(c) . If 'jn = 2 and holds, then fin {&,R) such that Z = iie*® has density 

2 ^\f*i\^\) < 1^1 < Further, let fin be as in ([13]) with a„ = (n^"i Vin = 2 

and ili. ji|) holds, then fi*^ fi* with fi* = Unif{\z\ = 1}. 

From the proof of Theorem [S] it is actually seen that the condition is differen¬ 

tiable” is not necessary in (a). The general conclusion is that fin fi where fi is the product 
measure of Unif[0, 2 tt] and the probability measure on [0, oo) with cumulative distribution 
function F~^{x). We write the way in (a) to avoid a lengthy statement. In particular, (a) 
is general enough to our applications next. 

The values of Uj’s in Theorem [3| can be very different. Now let us play them and find 
out their limiting distributions. The first one below is on the size n^’s that are at the same 
scale and m is fixed. 


Corollary l Let mn = m > l and lim^^oo ^ = Oij ^ [0, 1] for 1 < j < m. Assume 
Zi, • • • , Zn are the eigenvalues ofiK^'^\ Let a„ = (DI^i Ln be as in 

J I'j 

(1) . If ai = ■■■ = am = 1, then fin Lfnif{\z\ = 1}. 

(2) . If ai = ■ ■ ■ = am = a G [0,1), then fin ^ it* with density ^ \^\ ^ 

where f*{x) = 2(1 — for x G [0,1]. 

(3) . If aj < 1 for some I < j < m, then fi^ fi* with density 2 ;^/*(kl) on the set 
{0 < \z\ < 1} where f*{x) = ^F~^{x). 

Trivially, part (2) in the above corollary is a special case of part (3). We single it out since 
f*{x) has an explicit expression. Picking m = 1 in ( 1 ) of Corollary (T] we know that, with 
probability one, fi^ Umf{\z\ = 1}, where < = ^ ^Zj/a^ with a„ = (n/ni)^/2 i 
as n —)• oo. It implies that ^ converges weakly to Unif{\z\ = 1}. This conclusion 

is obtained by Dong et al. (2012). Taking m = 1 and ai G (0,1) in (2) of Corollary [3 
we get a result by Petz and Reffy (2005): ^Yl]=ibzj converges weakly to a probability 

measure with density f{z) = for \z\ < y/af (noticing the scaling in fi*n is different 

from 

Pick a = 0 from (2) of Corollary d] the limiting density becomes ^ for l^l < 1 , 

which is exactly the same as that of the product of Ginibre ensembles; see the paragraph 
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above Theorem [2j This is not a coincidence. In fact, Jiang (2009) show that the n x n 
submatrix Xi of the ni x ni matrix Ui can be approximated by a Ginibre ensemble as n = 
o{y/ni) in the variation norm. Similar conclusion also holds for Haar-invariant orthogonal 
matrices (Jiang, 2006). 

If m = nin oo, and {uj; 1 < j < m} are almost sitting on a curve, what is the 
corresponding limit appearing in Theorem [3l? To answer the question, assume there exists 
a continuous function q{x) defined over [ 0 , 1 ] satisfying 0 < q{x) < 1 for 0 < x < 1 and 



( 1 . 12 ) 


Define 


F{x) = xexp ( “ / log (l “ 9(^)(1 “ ® 


fix) 


mf ,t 

X Jo ^-Q{t)i^-x) 


for 0 < X < 1 and /(O) = F{0) = 0. 


(1.13) 

(1.14) 


Corollary 2 Assume m = mn —?> oo and IJ.iHI) holds. Let F{x) and f{x) he as 
in lil.lA) and 1)1.14^ , respectively. Let Zi,--- ,Zn be the eigenvalues of Set bn = 

YljLiin/nj) and 

1 


n 


i=i 


Then ^ (0, R) and Z = iZe*® has density 27 r|z| /(F-ThD) ® ^ ^ 


For an Haar-invariant unitary matrix, the empirical distribution of its eigenvalues is 
asymptotically the arc law Unif{\z\ = 1}, see, for example, Diaconis and Shahshahani 
(1994) and Diaconis and Evans (2001). If n is very close to nj for each j in Theorem [31 
that is, the truncated sub-matrix with size n x n of Uj is almost the same as Uj for each 
j, do we always expect the arc law Unif{\z\ = 1}? The answer is no and, as a matter of 
fact, it depends on the sum of Ij = nj — n for 1 < j < m. 


Corollary 3 Let m = mn —>■ oo as n ^ oo and lim^^oo maxi<,<m„ I— — 1| =0. Let 

—J — rij 

lim^^oo ^ = /3 G [0, oo]. Assume Zi,--- ,Zn are the eigenvalues ofLet /x* 

be as in [T^ with o„ = (Oj^i 

(a) . If 13 = 0, then ^ Unif{\z\ = 1}. 

Let Hn be as in (|1.3p with hn{r) = (r'^fbn)^^'^'^ and bn = 

J Vj 

(b) . If 13 G (0, oo), with 7 „ = 2 we have pinh with density ,r|^p(/ 3 -^iog|z|)^ f^'^ 0 < j^j < 1. 

(c) . If (3 = oo, with 7 n = 7 YlT=ih 'xe have pin T where p, has density 2 -K\z\^(v-\og\z\f 
for 0 < \z\ < 1. 
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Finally, we work on the case that n is much smaller than n^-’s. 

Corollary 4 Let m = mn —t OO and maxi<j<m„ ;^ = 0 as n —>■ oo. Review Zi , - ■ ■ ,Zn 
are the eigenvalues of . Set K = IljLiin/nj) and 

1 

i=i 

Then (©j-R) o-nd Z = y/Re^® follows Unif{\z\ < 1}, that is, the circular law. 

Picking m = 1, since Zj = |Zj|e*®-’, by the continuous mapping theorem, we get that, 
with probability one, \ 'YTi=\^{njlnYI'^Zj converges weakly to the circular law Unif{\z\ < 
1}. This result is found and proved by Dong et al. (2012). 

1.3 Structures of Determinantal Point Processes on Complex Plane 

In this section we state our results on rotation-invariant determinantal point processes on 
complex plane; see the set-up in Lemma ll.il below. The normalizing constant of their joint 
density function, moments and the structures of the two product matrices aforementioned 
are obtained. 

Let {Zi,--- , Zn} be n complex-valued random variables. Let K{z,w) : —)• C with 

K{z, w) = K{w, z) for all (z, w) € C^. Let be a Borel measure on C. We say {Zi, • • • , 
forms a determinantal point process with kernel K(z, w) and background measure n if the 
density function of {Zi, ■ ■ ■ , Z^j is given by 

fn(zi,--- ,Zk) := det{K{zj,zi))^^ {zi, ■ ■ ■ , Zk) e , (1.15) 

with respect to the product measure for all 1 < A; < n. The determinantal point process 
given above is a special case of a general dehnition in which the space C can be a arbitrary 
metric space. The dehnition here is good enough for our discussion. For general case, one 
can see, for example, Soshnikov (2000), Johansson (2005) or Hough et al. (2009) for a 
reference. 

Let <^(x) > 0 be a measurable function dehned on [0, oo) with 0 < dy < oo 

for each 1 < j < n. Dehne 

_ y2i-i(^(y)j(y > 0) 

fo"y^^-My)dy ^ ^ 

for 1 < j < re. Dehne 

1 ” 

Pn{y) =-'^Pj{y), y>0. (1.17) 

re 

i=i 

Now we start a series of results on the determinantal point processes. The following 
is a special case of Theorem 1.2 from Chafai and Peche (2014). It is another version of 
Theorem 4.7.1 from Hough et al. (2009). 
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Lemma l.l Let ^p{x) > 0 6e a measurable function defined on [0,oo). Let ,Zn) 

be the probability density function of (Zi,--- ,Zn) G C” such that it is proportional to 
ni<j<fc<n ki “ Zj\). Let Yi, - ■ ■ ,Yn be independent r.v. ’s such that the density 

of Yj is proportional to y'^^~^(p{y)L(y > 0) for each 1 < j < n. Then, g{\Zi\, • • • , \Zn\) and 
g(Yi, ■ ■ ■ ,Yn) have the same distribution for any symmetric function g{yi, ■ ■ ■ ,yn)- 

The next result is mostly known. Our contribution is that we are able to evaluate C for 
any ip{x), where f{zi, ■■■ ,Zn) = C ■ ni<j<fe<n ki “ is as in LemmalMl 

Proposition l Let ip and f be as in Lemma \1.1[ Set = 27r x^^~^^p(x) dx for 

all k = 0, ■■ ■ , n — 1. Then, C~^ = uIcqCi ■ ■ ■ c^-i and (Z^, ■ ■ ■ , Zn) forms a determinantal 
point process with background measure p{\z\)dz and kernel K{z,w) = 

The next result gives an estimate of the fourth moment of the sum of a function of Zj's, 
where Zj^s forms a determinantal point process. Hwang (1986) obtains a similar result for 
the special case of the complex Ginibre ensemble with p{z) = In particular, we do 

not assume any differentiability of p{x). 


Proposition 2 Let Zi, - ■ ■ , Zn and p{x) be as in Proposition{J\ with p{dz) = p{\z\)dz. 
Then, for any measurable function h{z) : C —)• M with sup^g^ \ ^ have 


E[Y^ih{Zj)-EhiZfi) 

i=i 


< Kn'^ 


for all n > 1, where K is a constant not depending on n, p{z) or h{z). 


The essential of the proof of Proposition [2] is the estimate of E\^-^fih{Zj) — Eh{Zj)). It is 
carried out by using ()1.15l) repeatedly. Our proof is different from the analysis of a contour 
integral by Hwang (1986), which seems to ht the Gaussian kernel p{z) = only. 


Proposition 3 Let Zi’s and p{x) be as in Lemma \l.l[ Then the following hold. 

(i) Let Pn{ ) be as in d 1. j7| ). Then, for any bounded measurable function h{z), 

Eh{Zi) = ^ (^l\{rP<^)de)Pn(.r)dr. 

(a) Let 01 be as in Then, \Zi\ has density Pn{r), 0i ~ C/m/[0, 27r] and the two are 

independent. Consequently, for any bounded measurable function g{r,9), 

E5(0i,|Zi|) = ^ (^J^\{9,r)d9)Pr,ir)dr. (1.18) 

The following result reveals the structure of the eigenvalues of the product of Ginibre 
ensembles. The key is the Gamma distribution. We will switch the use of the notation of 
“r” in the next two lemmas, which will serve as an index instead of the radius of a complex 
number used earlier. 


9 








Lemma 1.2 Let Zi, - ■ ■ , have density as in Let {sj,r) l<r< m} 

be independent r.v. ’s and sj^r have the Gamma density y^~^e~yI{y > 0)/{j — 1)! for each j 
andr. Then g{\Zi\‘^, ■ ■ ■ ,|^nP) and giUfL^ si^r, ■ ■ ■ ,11" Li Sn,r) have the same distribution 
for any symmetric function g{ti, ■ ■ ■ ,tn)- 

Recall the beta function 



(1.19) 


a > 0, 6 > 0. 


The following lemma describes the structure of the eigenvalues of the product of trun¬ 


cations of Haar-invariant unitary matrices. It has the same setting as Lemma 11.21 with 


“Gamma distribution” replaced by “Beta distribution”. 

Lemma 1.3 Let Zi, - ■ ■ , Zn have density as in Q. Let {sj^r, l<r< m} 

be independent r.v.’s and sj^r have the Beta density < y < 1) 

for each j andr. Then g{\Zi\‘^ ,■ ■ ■ ,|Z„p) and g{Y\^^i si^r, ■ ■ ■ ■,W!^=iSn,r) have the same 
distribution for any symmetric function g{ti, ■ ■ ■ ,tn). 

Comments. We now present some remarks and state problems for future. 

1. There are other type of studies on the product of random matrices in literature. 
The size of each matrix is assumed to be fixed and conclusions are obtained by letting the 
number of matrices go to infinity. Two typical interests of the product matrices are their 
norms and entries; see, for example, Furstenberg and Kesten (1960) and Mukherjea (2000). 

2. In this paper we study two kinds of product matrices: the product of Ginibre 

ensembles and that of truncated Haar unitary matrices. Notice the Ginibre ensemble and 
truncated Haar unitary matrices are of the Haar-invariant property. We believe that the 
same method (Theorem [T|) can be used to derive the spectral limits of the products of 
other independent Haar-invariant matrices. The key is the explicit formula of ip{x) and the 
verification of ‘Vn as stated in Theorem [TJ 

3. The universality of Theorem [2] is an interesting problem. Namely, replacing the 
normal entries in the Ginibre ensemble with i.i.d. non-Gaussian random variables, does 
Theorem [2] still hold? In fact, Bordenave (2011), O’Rourke and Soshnikov (2011) and 
O’Rourke et al. (2014) show that this is true for fixed m. We expect Theorem [2] to hold 
for non-Gaussian entries and arbitrary m = mn- 

4. There are three Haar-invariant matrices generating the Haar measure of the classical 
compact groups: Haar-invariant orthogonal, unitary and smplectic matrices; see, for exam¬ 
ple, Jiang (2009, 2010). Similar to Theorem [3] one can work on the same limiting problems 
for the orthogonal and symplectic matrices. 

5. If we change the square matrices Xj in Theorem [2] to rectangular matrices and keep 
the Gaussian entries of each matrix, that is, Xj is Uj x Uj+i with Um+i = ni, it will be 


10 





interesting to see the corresponding result. The limiting distribution will have a rich feature 
as the ratio nj/rej+i fluctuates for each j. 

6. Let Zi^ - ■ ■ ,Zn be the eigenvalues of the product OJli where Xi, • • • , are 
i.i.d. Ginibre ensembles. The spectral distribution of Zi, • • • , is well understood through 
Theorem [21 The transitional phenomenon of the spectral radius maxi<j<m \ Zj\ is obtained 
by Jiang and Qi (2015). It is classified by c := lim ^ with c = 0, c G (0,oo) and c = oo. 
The spectral radius of X^™) in Theorem [3] with m = 1 is also investigated in the same 
paper. With the help of Lemma 11.31 the spectral radius of X*-™^ for arbitrary m can be 
done similarly. 

The rest of the paper is organized as follows. In Section 12.11 we prove the results 
stated in Section [1.31 that serve as technical tools to prove the main results. We then prove 
Theorems [T]l3] and all of the corollaries in Sections 12.2112.41 


2 Proofs 

In this section we will prove the main results stated in the Introduction. We first prove 
those in Section oi since they serve as tools to derive the main limit theorems. Theorems 
m and all corollaries are proved one by one afterwards. 


2.1 Proofs of Propositions [I}[3], Lemmas 11.21 and 11.31 


Proof of Proposition [T1 For 1 < j < n, let Zj = with Rj > 0 and Qj G 

[0,27r). The Jacobian is obviously equal to Thus the joint density function of 

,Rn,Qi,--- ,&n) is given by 


n 


( 2 . 1 ) 


Set 


AR = 


( 


T\e 


i6\ 


1 


r2e 


i02 


1 \ 


1)02 . . . j.n—l^i{n—l)0n^ 

By the formula of the Vandermonde determinant, the hrst product in (12.ip is equal to 
(det(M„))^. It follows that 


where 


Yl - rfce *'^'=|2 = C/V 

l<j</c<n 


U = det(M„) = ^ sign((T) 

u j = l 
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and cr is a permutation running over all elements in the symmetric group Sn- 

Note that Jjq = 0 for any integer j ^ 0. Therefore, any two terms in the sum 

are orthogonal to each other, that is, for any a ^ o', 


Thus, 



i=i 



UUdOi 


n 

den = {2T^TY.X\^T^'^~^^- 

a j = l 


By integrating out all dj’s in (12.1|) . we get that the probability density function of (Ri, • • • , Rn) 
is equal to 


n n 

C ■ (2,r)”(ri ■ ■ ■ r„) n = C ■ (2^)” J] 

O-eSn i = l 0-£Sn j = l 

for ri > 0, • • • , > 0 and the density is 0, otherwise. It follows that 


1 = 


~ n 

c(2^r / E n 

n 




^p{rj) dr I ■ ■ ■ dr-n 


c{2^r E / n 


^(f{rj)dri-■ ■ drn- 


For each cr S S'n, it is easy to see that the integral is equal to ^(p{x) dx. We 

then get the value of C. 

Second, the first step says that the density function of {Zi, • • • , Z„) is 


,Zn) = ^ 


n\ Cq ■ ■ ■ Cn—l 




Now write 


(co • • • Cn-i)i/2 


{zj - Zk) = det 


l<j<k<n 


l<ji<fc<n 

J=1 


^ Poi^i) 

P0{Z2) 

Po{Zn) ^ 

Pl{zi) 

Pliz2) 

PliZn) 


Pn-l{z2) •• 

Pn—l{^n) j 


where pi{z) = z^ jyfcl for I = 0,1, 2 • • •. Let A be the above matrix. Then, 
1 


^0 * * * C 77 ,—1 


JJ \zj — Zk\‘^ = det(A*A) 


l<J<fc<Tl 


n—1 


( 2 . 2 ) 


. ^ Pkj^ijpkjz^j)') — det{[K (^Zi, Zj)nxn) ■ (2-3) 


k=0 
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By the polar transformation, f ^2 z^(p{\z\) dz = 27r dr = Cj for j = fe, and 

the integral is equal to 0 for any non-negative integers j / k. Hence {po{z),pi{z),p 2 {z), ■ ■ ■} 
are orthonormal with respect to the measure dz. By using Exercise 4.1.1 from Hough 

et al. (2009), we get the desired conclusion from (j2.2p and (j2.3p . ■ 

To prove Proposition [21 we need some basic facts regarding point processes. Let 
Zi,--- ,Zn be random variables with symmetric density f{zi,--- ,zn) with respect to 
reference measure p on C. View them as a point process ^Zi- Then this process has 
n-point correlation function 


Pn{,Z \) ■ ■ ■ ) Zri) — 


N\ 


[N — n)\ 


f{zi, ■■■ ,Zn) p{dZn+l) ■ ■ ■ p{dZN) 


for 1 < n < V and p]\f{zi,--- ,zn) = N^-f{zi,--- ,zn)- Let fnizi,-'' be the joint 
density of Zi, ■■ ■ , Zn- Then, 


* * * ; ^n) — jyi Pny^l-i * * * •, 


(2.4) 


for 1 < n < V. This means that, for any measurable function h{zi, • • • , Zn) with 1 < n < V, 
we have 

i N f ^ 

Eh(^Zi, • • • , Zn} — I h(^zi, ■ ■ ■ , Zfi)Pn{zi 1''' 1 Zn) I I pidzj) (2-5) 

IV! Jc" 


provided 


„ n 

/ \h{zi,--- ,Zn)\Pn{zi,--- ,2n)]T p{dZj) < OO. 

Jc^ ,=1 


See further details from, for example, Johansson (2005) and (1.2.9) from Hough et al. 
(2009). Let Zi, - ■ ■ , be as in Proposition [H Then 


n—1 


Pkizi,-- ■ ,Zk) = det {K{zi,Zj))^^. _ where K{z,w) = ^ —{zw) 


( 2 . 6 ) 


k=0 


Evidently, K{z,w) = K{w,z) for all {z,w) G C^. Eurther, the product of the diagonal 
entries of {K{zi, Zj))i<ij<k is equal to 0^=1 Pi(^l)- 

Proof of Proposition [2l Eor convenience we now switch the notation n to N. So we 
need to prove 

1 4 


E 


N 


< kn"^ 


(2.7) 


i=i 


for all V > 1, where iL is a constant not depending on N, ^{z) or h{z). 
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Obviously, h{z) = {z)—h~{z) where:= max{/i(z), 0} and/i“(z) := min{/i(z), 0}. 

With the trivial bound (o + < 8(a^ + 6^) for (a, 6) G to prove the proposition, we 

can simply assume 0 < h{z) < 1 for all z G C. 

Set eh = E/i(Zi). Then en G [0,1]. Without loss of generality, we assume /x is a 
probability measure, i.e., /x(C) = 1. Review the identity that 

(ai + • • • + 

N 

= + 12 ^ a^ajafc + 24 ^ aiUjUkai 

1=1 i<j i^j i^j^k i<j<k<l 

where are all between 1 to N. Set aj = h{Zj) — e/^, 1 < j < N. Since h{z) is a 

bounded function, it is easy to see that, to complete the proof of (I2.7p . it suffices to show 

2N{N - 1){N - 2)E[(h(Zi) - eh)\h{Z2) - eh){h{Zs) - e^j < OiV^ ( 2 . 8 ) 


and 


N{N - 1)(1V - 2){N - 3)E[l[{h{Zj) - Bh) 

1=1 


< 82iV^ + 18iV 


(2.9) 


for all iV > 4. Define Uk{zi,--- , Zk) = 0^=1 Pi(^l) “ Pk{zi., ■■■ ,Zk) for any A: > 1. By using 
the Hadamard inequality, Pk{zi, • • • , Zk) < 0^=1 know Uk{zi, • • • , Zk) > 0. Then 

using (12.4p 


Ukizi,--- ,Zk) 


N' 




1=1 


N{N-l)---{N-k + l)fk{zi,--- ,Zk). 


Since 1 — nj=i(l ~ ^j) ^ Sj=i S [0) 1]) 1^6 obtain 

f Ukizu ■■■,Zk)T\ p{dzj) =N’^-N{N-l)---{N-k + l)< N’^-K (2.10) 

It follows from (|2.4I) that 


2NiN - 1)(1V - 2)E[(h(Zi) - eh)\h{Z2) - eh){h{Z^) - e^j 
2 / {h{zi) - eh)'^{h{z 2 ) - eh){h{z3) - eh)p3izi,Z2,Z3)Y\_p{dzj) 


1=1 


= 2 


f ^ 

/ {h{zi) - ehf{h{z 2 ) - eh){h{z 3 ) - eh) pi{zi) pi{z 2 ) pi{z 3 ) TT p,{dz 
Jc^ 

r ^ 

2 / {h{zi) - eh)^{h{z 2 ) - eh){h{z 3 ) - 6/1)^3(21,22,23) TT p.{dzj). 

f=i 
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Since pi{zi)pi{z 2 )pi{z 2 ) = h{zi)h{z 2 )fi{z 2 ,) from ([23D, we have 

2 / {h{zi) - ehf{h{z 2 ) - eh){h{z 3 ) - eh)pi{zi)pi{z 2 )pi{z 3 ) TT p{dzj) = 0. 

Thus, we get from (|2.10p and the assumption h{z) € [0,1] for all z G C that 
2N{N - 1){N - 2)E[(h(Zi) - ehf{h{Z2) - eh){h{Z3) - e0]| 

f ^ 

-2 / {h{zi) - ehf{h{z 2 ) - eh){h{z3) - eh)u{zi,Z2,Z3)Y\_fJ-{dzj) 


/C3 


j=i 


< 


f ^ 

2/ u{zi,Z 2,Z3)Y\ p{dZj) 

Jca fJi 


= 6N^, 


proving 


Now, we start to prove ()2.9p . 


Define = E[n ■ h{Zj)] for 1 < /c < 4. Note that 


ll{h{Zj)-eh) = et-elY. h{Z,) + el hiZi)h{Z,) 

j=l l<i <4 l<*<i <4 

-eh h{Zi)h{Zj)h{Zk) + h{Zi)h{Z2)h{Z3)h{Zi). 

l< 2 <_ 7 <fc <4 

Taking expectations on both sides and noting that bi = eh we have 

4 


E[n(h(Z,-)-e,) 


= -Set + 6e?62 - delfts + 64. 


i=i 


It follows from (12.51) that for 2 < /c < 4 


h = 


{N -k)\ 

m 

{N -k)\ 

m 

{N -k)\ 


K K 

,Zk)'[[p{dzj) 

J=1 

k k 

YlKdz 


i=i 


C'' S=1 


i=i 


+ 


m 


/c'= 




1=1 


zir-- ,Zk) -Wpi{zj)) Wp{dzj) 

j=i j=i 


{N-k)\N'^ ^ {N-k)\ 

'^h ~ 


N\ 


N\ 


[flhiz 




Uk{zi,--- ,Zk) Y\_p{dzj 
i=i 


We next evaluate the last integral for k = 2,3, 4. Denote 

k k 


“fc = / TT dizj) Ukizi, • • • , Zfc) TT pidzj). 


( 2 . 11 ) 


( 2 . 12 ) 
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Then we have 


(2.13) 


, (N-kV.N'^ . (N-k)l 

h = - -TTT- 4 - ak 


for k = 2,3, 4. Obviously, > 0 for 2 < /c < 4 since Uk{zi, ■ ■ ■ , Zk) > 0 for all (zi, - ■ ■ , Zk) G 
C^. In order to prove (I2.9jl . we need to study b 2 ,bs, 64 in (I2.11h . Based on (I2.13|) . it suffices 
to work on a 2 , as, a 4 . We will do so step by step in the following. 

Estimate 0 / 0 : 2 . It follows from (I2.10h that 


(2.14) 


0 < a2 < N. 


Meanwhile, by the determinantal formula (j2.6l) . since that 


P2{zi,Z 2) = K{zi,Zi)K{z2,Z2) - K{zi,Z2)K{z2,Zi) 

= Plizi)pi{z2) -\K{zi,Z2)\‘^, 


we have 


/ h{zi)h{z2)\K{zi, Z2)\‘^p{dzi)p(dz2) 
Jc^ 


(2.15) 


02 = 


which will be used later. 


Estimate of as- Set 



(2.16) 


In this step, we will show 


maxjoa, l/^a]} < 3IV^. 


It is easily seen from ()2.10l) that 03 < 3iV^. We now estimate f3s- By (12.61) again, 

us{zi,Z2,zs) 

= pi{zi)pi{z2)piizs) - Psizi,Z2,Zs) 

= K{zi,Zi)K{z 2 ,Zs)K{zs,Z 2 ) + K{zi,Zs)K{z 2 ,Z 2 )K{zs,Zi) + K{zi,Z 2 )K{z 2 ,Zi)K{zs,Zs) 
-K{zi, Z 2 )K{z 2 , Zs)K{zs, Zi) - K{zi,Zs)K{z 2 ,Zi)K{zs,Z 2 ) 


Pl{zi)\K{z 2 , Zs)\^ + pi{z 2 )\K{zi,zs)\‘^ + pi{z3)\K{zi,Z2)f 
-K{zi, Z 2 )K{z 2 , Zs)K{zs, Zi) - K{zi,Zs)K{z 2 ,Zi)K{zs,Z 2 ). 


(2.17) 


All three functions in (I2.17j) are nonnegative. For the first term in (I2.17j) we have from 


16 

















(|2.15p and then f|2.4p that 


« c$ CS 

/ Pl{zi)\K{z2,Z3)\‘^Yln{dZj) 

i=i j=i 

Jc JC2 

02 / /i(zi)pi(zi)/i(d2:i) 

Jc 


= Na 2 eh. 

The same is true for other two terms in (j2.17l) . Trivially, 


/C3 


/^3 / tt ^ ^ ^ ^ 7 ^3 ) (-^2 7 ( Z3 , -^2 ) n lJL{dZj). 

Therefore, we obtain 


i=i 


03 = 


/C3 


o o 

[ n K ^ j ) ^3(21, 2:27 23) Jt p { dzj ) = 2>Na2eh - 2 / 33 , ( 2 . 18 ) 


i=i i=i 

which together with the facts a 2 < N and 03 < 2>N‘^ implies 

lAI < f +< 3iv^. 


Estimate of a^. This step is a bit involved. The sketch of the proof is as follows. Since 
P 4 (zi, 227 23 , 24 ) is the determinant of {K{zi, Zj)) - it can be written as the sum of 24 
terms: 


P4{zi,Z2,Z3,Z4) ='^ sign{a)K{zi, Z„(^ 1 ))K{Z 2 , Z„(2))K{z3, Za(3))K{z4 , 2 ^( 4 )) 7 

a 

where a = ((t(1), (t(2), (t(3), (t( 4)) runs over all 24 permutations of (1,2, 3,4). Excluding 
(1,2,3,4), all other 23 permutations can be classified into one of the following 4 sets: 

Di = {(1,2,4,3), (1,4,3, 2), (1,3,2,4), (4,2,3,1), (3, 2,1,4), (2,1, 3,4)}, 

3^2 = {(2,1,4,3), (3,4,1, 2), (4, 3, 2,1)}, 

D 3 = 1(4,1,2, 3), (4, 3,1, 2), (3,4, 2,1), (3,1,4, 2), (2, 3,4,1), (2,4,1, 3)} 

and 


A = {(I 7 4, 2,3), (1,3, 4, 2), (4,2,1, 3), (3, 2,4,1), (4,1,3,2), (2,4,3,1), (3,1, 2,4), (2, 3,1,4)}. 
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Define 


Tk{zi,Z 2 ,Zs,Zi) = K{Zi^Mi))K{z2,Z^(^2))K{z3,Z„(^3,))K{z4,,Z^(^4^)) 

aeDk 

for k = 1,2,3,4. Then 

4 

P4{ZI,Z2.,Z3,Z4) = K{zi,Zi)K{z2,Z2)K{z3,Z3)K{za,Z 4) + '^{-l)^Tkizi, Z 2 , Z 3 , Z 4 ,). 

This implies 


k=l 


Ui 


{zi,Z2.,Z3,Zi) = y~](-l)^+^rfc(zi,Z2,^:3,^4). 


(2.19) 


k=l 


Within each class Dk, all K{zi, Zf^(^i))K{z2, Zf^(^2))I'^{zs, Zf^(^3))K{z4,, Zf^i^4)) contribute equally 
to the integral [ 0^=1 H^j)] Tk{zi, Z 2 , Z 3 , Z 4 ) 0^=1 p{dzj). We have 

4 4 

[ n Kzj)]Ti{zi,Z2,Z3, Z4) JJ pidzj) 


>C^ ,=i 


1=1 


4 4 

6 / [Y\h{Zj)]K{zi,Zi)K{z2,Z2)K{z3,Z4)K{z4,Z3)Y\p{dZj) 

1=1 1=1 

4 4 

® L, ^3)Y[Kdzj) 

1=1 


,=1 


= P>N‘^e\a2 

by using (I2.15P and the definition that Ch = E/i(Zi). For r 2 , we have 
^ 4 4 


[ n ^2(^l> ^2, Z3, Z 4 ) JJ p{dZj) 


,=1 


1=1 


3 / [ J^/i(2:j)]ii:(zi,2:2)i^(^2,2l)i^(2:3,2:4)i^(2:4,^3) 


'C4 ,=1 


1=1 


= 3 


/ h(zi)h(z2)|i^(2:i, 2;2)Pli(d2;i)/i(d22) 

Jc^ 

/ /i(2;3)h(z4)|ii:(z3,2:4)p^(d2;3)/i(dz4) 
JC2 


/C2 

X 


= 3a^ 
< SiV^ 


by (I2.14p . Noting that 0^=1 i® ^ probability measure, we have from the Cauchy- 
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Schwarz inequality and the fact 0 < h{z) < 1 for all 2 ; G C that 


/c-i 


[ n ^( 2 j)]T 3 ( 2 i, 22,23,2:4) 


i=i f=i 

4 4 

6 I [Y]_h{Zj)]K{zi,Z4)K{z2,Zi)K{z3,Z2)K{z4,Z3)Y]_n{dZj) 

f=l 
1/2 


'c^ /=1 


< 


„ 4 4 

i=i i=i 

/ /i(2;i)/i(2;4)|ii:(2;i, 2;4)l^//((i2:i)/i(dz4) 

JC2 J 


1/2 


= 6 
= 6a^ 

< eiv^ 


by (I2T^ and (imH . 

We next estimate the term on r 4 . In fact, 


/C'l 


[ n H^j)] Tiizi, Z 2 , Z3, Z 4 ) fJ-idZj) 


i=i 


/=i 


8 / [ JJ/i(2;j)]K(2;i,2:i)K(2:2,Z4)K(2;3,2:2)K(2:4,2;3) JJ/i(d2;j) 


'C^ ,=i 


i=i 


4 4 

= 8/ [ JJ/i(2:j)]/9i(2;i)ii:(2;2,2;4)ii:(2;3,2;2)-fC(2;4,2;3) ]^/i(d2;j) 

i=i i=i 

4 4 

= SiVCft f ['^h{Zj)]K{z2,Z4)K{z3,Z2)K{z4,Z3)Y^I^{dZj) 
j=2 j=2 

= 8Nehl33, 


by due]). 

Now multiplying 0/=! h{zj) on both sides of (j2.19l) and integrating with respect to the 
measure n,=i IJ^idz^ ,) we obtain from (12.121) that 

Q!4 = QN'^e\a2 — 8Nehf33 + d23, (2.20) 
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where \d 23 \ < 9N‘^. From (I2.13p . (I2.18P and (12.201) we see that 

. ^ _ «2, . 

^ A^-1 iV(iV-l)’ 

, _ N'^el _ 3Neha2 - 2/33 

^ “ (A^ - 1)(A^ - 2) N{N - l)(iV - 2) ’ 

_ N^el QN‘^e\a2 - SNehPs + ^23 

^ “ (AT - l)(Ar - 2)(iV - 3) N{N -1){N - 2){N - 3) ’ 

Then it follows from ^2.111) that 


E[ll{h{Z,) - en) 

i=i 

= -3e^ + 6e|62 - de/^fts + 64 

_ / ^ 6A^ 4A^2 ^3 . 

V A^ - 1 ~ (A^ - 1){N -2)^ {N - 1){N - 2){N - 3)/ 

/ 6 12A^ 6A^2 . ^ 

V N{N - 1) A^(A^- l)(A^-2) “ N{N - 1)(A^ - 2){N - 3)/^'^“^ 

_ ( _?_^^ e.B. 

VA^(A^ - 1)(A^ - 2) A^(Ar-l)(A^-2)(A^-3)/ ^ 

d23 

~ N{N - 1){N - 2){N - 3) 

_ {3N + 18)e^ (6A^ + 36)e^a2 

(A^- l)(A^-2)(A^-3) ~ Ar(A^-l)(Ai-2)(A^-3) 

_ 24efe/33 _^23_ 

A^(A^ - 1)(A^ - 2){N - 3) N{N - 1)(A^ - 2){N - 3) 

90A^2 54JY 

- A^(Af- l)(A^-2)(A^-3) 

for A^ > 4 by the facts 0 < < 1, 0 < 02 < |/33| < 3N‘^ and |d 23 | < 9A^^. This proves 

P2.9p . The proof is then completed. ■ 


Proof of Proposition\^ We will only need to prove (ii). In fact, conclusion (i) follows from 
(ii) since h{z) = h{\z\e'’^) with 2 : = | 2 ;|e*®. By Proposition [D and p2.4l) . the density function 
/i( 2 i) of Zi is given by 

h(zM\A) = -K(z, ?Mki) = IE MihM. 

n n cu 

k=0 

Write z = X + yi = re*® with r > 0 and 6 G [0,27r). Then {x,y) = (rcos0,rsin0). The 
Jacobian for the transformation is known to be r. By (|1.16l) . (|1.17p and Proposition [TJ this 
implies that the joint density function of |Zi| and 0i is given by 


1 ^ r2‘^(r) 1 r2‘+V(r) 1 „ , , 

= — = - 2 .-:- = l^Pn(r) 


fc =0 


Cfc 


k=0 


C-k 


27r 
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for r > 0 and 6 G [0,27r). Therefore, \Zi\ and 0i are independent, the density function of 
\Zi\ is Pn{r), and 0i is uniformly distributed over [0,27r). The conclusion (|1.18l) follows 
immediately. ■ 

Proof of Lemma 11.21 Let Ynj, 1 < i < u, be independent random variables such that 
Ynj has a density function proportional to for each j. Then it follows from 

LemmaOthat , l^nP) and g{Y^i, ■ ■ ■ ,Y^J are identically distributed. In the 

proof of their Lemma 2.4, Jiang and Qi (2015) show that Y^j has the same distribution as 
that of ^j,r for each 1 < j < n. This yields the desired conclusion. ■ 

Proof of Lemma 11.31 Define = (1 — sY ^/(O < s < 1) for I > 1. Then we have 

Wi\s) = {I/T^). Set Pi^\s) = p^’‘^\s) and define Pm' recursively by 

.s^du 
u 


Jo ^ 


( 2 . 21 ) 


for m > 2 and positive integers h, - ■ ■ , Im- Evidently the support of pm’ is [0,1]. 

By induction, it is easy to verify from (jl.Sp that 


- m 

TT 


( 2 . 22 ) 


i=i 


for m > 1. Dehne 


0ai,-,M(i) = f gt ^p<Y’-M{s)ds, t>0. 

Jo 

0bviously, Now, for any m > 1, we have from (I2.2ip that 


1 rl 


) Jo u u 

Vm-I du. 

Keeping in mind that the support of pm’ ’*"*^( 5 ) is [0,1], the above is identical to 

I' v<‘”>{«) ([ s'-'pSir'”-*’ 


/ 

f 


_ ad 

- 


Jo 


'0 
. 4-1 




1 


du 




We thus conclude from the recursive formula that for any m>\, 


rn 

do r=l 


(2.23) 
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Let qj{y) be a density function proportional to for 1 < j < n. 

Thus, qj{y) is also proportional to y^~^Pm' '^"'\y) from ()2.22p . that is, for some Cj > 0, 
qj{y) = c~^y^~^pm’ Let Ynj be a random variable such that the density function 


of YY is qj{y). Since qj{y) is a density, then cj = y^ ^p^m’"'’‘""\y)dy 
nZimi r) from (I2.23p . and hence 


= a 


(^1 }■■■ >^m) 


Qj(y} = 


(h, 

y^ Pm 


0 <y <1. 


YC=lBU,lr) 

Denote by Mj{t) the moment generating function of logl^^-. Then, 

= / y^Qj{y)dy 

Jo 

foy^^*-^Pm’"'’'"‘^(y)dy 


Mj{t) = E(e 




n 


B{j + t, Ir) 


=1 B{j,lr) 


(J) = 


(2.24) 


for t > —j, which is the same as -^6xp{t log where {sj^r^ ^Yj<n, l<r< m} 
are independent random variables and Sj,, has the Beta density — y)^’'“^/(0 < 

y < 1) for each j and r. Then, Y^j and O^i ^j,r have the same distribution. This and 
Lemma o lead to the desired conclusion. ■ 


2.2 Proof of Theorem [T] 

Proof of Theorem [H Write Zj = with Qj S [0,27r) for 1 < j < n. We need to 

show that for any continuous function u{9, r) with 0 < u{9, r) < 1 for all 9 G [0, 27r) and 
r > 0, 

^ roo r 27 V 

-'^u{Gk,hn{Rk))^— / u{9,r)d9u{dr) a.s. (2.25) 

n ^ 27r 7o Jo 

as n —>• oo. Obviously, (0fc,iifc) for 1 < A: < n have the same distribution. First, by the 
Markov inequality and Proposition [21 

1 ^ 

P{\-'^u{Qk,hn{Rk)) - Euiei,hn{Ri))\ > e) 

^ k=l 

^ E[ELi {u{ek,hn{Rk)) -Eu{Gl,hn{Rl)))]^ 

C 

for every e > 0, where O > 0 is a constant not depending on n or e. This implies that 

OO 1 ^ 

{Gk,u{hn{Rk)) - E'u(0i, hn{Ri))) I > < oo. 

n=l k=l 
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We conclude from the Borel-Cantelli lemma that, with probability one, 


1 

- {u{&k, hn{Rk)) - Ett(0i, /in(-Ri))) -)■ 0 (2.26) 

^ k=l 

as n ^ oo. 

Note that G(r) = ^ u{6,r)d6 is bounded and continuous in r S [0, oo). Since Vn 

converges weakly to ly with probability one, we have 

1 ^ roo 

-y^G{hn{Yk)) ^ G{r)iy{dr) 

^ k=i 

with probability one. This implies 



I n 


n fc 

Y^G{K{Yk))- / 

fe=i 


G{r)v{dr] 


0 


via the bounded convergence theorem. Hence 


1 _ ^ _ /*oo 

-^y^G{K{Yk)) - G{r)iy{dr) ^ 0, 

^ k=i 

which together with Proposition [3] yields 


poo -I ^ poo 

Yu{Qi,K{Ri)) = G{hn{r))Pn{r)dr =-Yy^G{K{Yk)) ^ G{r)iy{dr 
Jo n Jo 


This and (|2.26l) imply (I2.25p . The proof is then completed. 

Now we prove the conclusion for (/r* , z/*,/r*, z^*). It suffices to show that, for any con¬ 
tinuous f{z) with 0 < f{z) < 1 for every z £ C, 


1 

n 




f{re''^)d0v*{dr) 


(2.27) 


with probability one. 

Define g{6,r) = f{rp^) and hn{r) = r/a^. By TheoremlH ()2.25p holds. It follows that, 
with probability one. 


1 

n 



completing the proof of (I2.27p . 


1 " 

^ ^ ^(Qfc) {Rk ) ) 

n 

k=l 

1 poo p2'K 

^ Jo Jo 

1 poo p2'K 

^ Jo Jo 
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2.3 Proof of Theorem [2] 

We first need a technical lemma as follows. 

Lemma 2.1 Suppose {/i„(x); n > 1} are measurable funetions defined on [0,oo) and Vn 
are defined as in do]). Let Yi, ■ ■ ■ ,Yn be as in Lemma [7771 and u be a probability measure 
on M. Then v if and only if 

n 

hm -y^P{hn{Yfi<r)=G{r) 

n—)-oo 77, 

for every continuous point r of G{r), where G{r) := z^((—oo,r]), r G M. 

Proof. Let Gq denote the set of all continuity points of G. Note that Un converges weakly 
to V with probability one if and only if r'„((—oo,r]) —>• z^((—oo,r]) with probability one for 
any r G Gq-, that is, for all r G Gq 

n 

-y^I{hn{Yfi<r)^G{r) (2.28) 

i=i 

with probability one. Since Li, • • • , W are independent random variables, 

1 ” 

z/n((-oo,r]) = - ^I{hn{Yj) < r), 

which is the average of n independent bounded random variables. By calculating the fourth 
moment, applying the Chebyshev inequality and then the Borel-Canteli lemma we can show 
that for any r G M, 

^ n 1 ^ 

- E ^ E s -•) 

i=i i=i 

1 " 

i=i 

with probability one. This and (I2.28jl imply the desired conclusion. ■ 

Proof of Theorem [2l Let hn[y) = y > 0. By applying Theorem [T] and 

Lemma O it suffices to show that 

n 

-^P(|Z,f/-Vn<y)^y, ye (0,1) 

which is equivalent to 

^ n mn 

<y)^y, y e (0,1) (2.29) 

j=l r=l 
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by Lemma 11.21 where sj^s are as in Lemma 11.21 Define r/(x) = x — 1 — logx for x > 0. 
Since t]{x) = it is easy to verify that 


0 < i](x) < 


2 min(x, 1) ’ 


X > 0. 


(2.30) 


Set Wj = n^i ^j,r for 1 < i < rr. Then, 


m„ 

log Wj = ^ log Sj^r 

r=l 


for each j. By using the expression logx = x — 1 — r]{x) we can rewrite log Wj as 

rrin ^ 

log Wj = y] log ^ + TTln log j 

r=l ^ 

^ rrin rrin 

= - -j)+"irilogi - 

^ r=l r=l 

Write Tj = for each 1 < j < n. Then Tj is the sum of jrrin i.i.d. random variables 

with the exponential distribution of mean 1. Hence E(rj) = and Var(Tj) = for 
1 < j < ri. From the above equations we have 



log Wj 


-{Tj - jiTin) + mnlogj - — 


r=l 


(2.31) 


Since Sj^r has the Gamma density ^I{y > 0)/(j — 1)!, the moment generating func¬ 

tions of log Sj^r is 

mj{t) = E(e*'°*^^-) = r y^y^-^e-ydy = 

J- Ul JO 


m 


for t > —j. Therefore, 




(2.32) 


The function ip is the so-called Digamma function in the literature. By Formulas 6.3.18 
from Abramowitz and Stegun (1972), 


1 


ipix) = logx — -—asx—)• -|-oo. 


2 x 


1 


Because E(loglTj) = hi log(sj,r) = ’mn'ipij), we see that 

m„ 

= mnilogj - tpij)). 


(2.33) 


(2.34) 


r=l 
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Now we fix y € (0,1). Write Sj = YlT=i vi^)- It follows from (j2.3ip that 


lo ^ II IIJI 

j=l r=l 

n 

P ^ log Wj < m, 

i=i 

n ^ 

'^p(^-{Tj -jrrin) + mnlogj - Sj < 


rUr 


i=i 


< y/j'^n log( —) + \ —Sj) 
^ ^y/jrrin J V V 


(2.35) 


For any fixed small number e G (0,1/2) such that y(l + e) < 1, define integers = 
[ny(l + e)] + 1 and j~ = [ny/(l + e)], where [x] denotes the integer part of x. Obviously 
we have 

ifj<Jn (2.36) 


nn 

^ > 1 + e 


and 


J 

ny 


< 


1 


if J > Jn- 


(2.37) 


j 1 + e 

Since Tj is the sum of jrrin i.i.d. random variables with both mean and variance equal to 
1. Then Var( (Tj — jmn)) = 1- From (|2.36l) and the Chebyshev inequality, 


/jnir, 


n 


Jn -| \ 1 1 


j P jmn(log(l+ e))" 

= o(l25!^)^o 

nrrin 


(2.38) 


as n —>• oo. This implies 


1 Jn 1 

lim - ^p(—:^{Tj - jrrin) < log (—)) = — 

1^00 n ^ yJimr, 7 / 1 


Note that Sj > 0 from ()2.30l) . We obtain by (I2.35P 


.. n rrin 

lim inf -V P((V\ Sjr)^/^^/n<y 

n—>oo fi V 


j=l r=l 
1 


1 * ^ 1 

> liminf-y~]pf — -.{Tj - jmn) < juin^ogj—)) 

n^Qo n \y/jmr, ^ 7 / 


1 + £ 


(2.39) 
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From (I2.37P we have 


^ E '='(7lb- v7^‘“i5(7) + J ^Sj) 


J=Jn 

n 




< 


n 




J=Jr7 


+ P(5,>^log(l+£)) 


1 n . - 1 

- n'^ - -2\^^^°s(l + e)) +- ^’(<5i > ^log(l + e)). 


J=Jn 


J —jn 


By the same argument as in (I2.38p . the first sum above goes to zero. Further, by the 
Markov inequality, p2.34p and then (I2.33p . the last sum is controlled by 


n ^ ^ log(l + e) 
i=iit 


^ j- y^ "in(log j - V^(j)) 
— ^ rn^ log(l + e) 

3=3t 

1 ” 1 

= 0(1) j: 1^0 

j=ji 


as n —)• oo since j = 0(1)- These and (I2.35h imply 


. n rrin 

limsup — ^ T’('( JJ jn < y 

n—>oo n _ \ . 


j=l r=l 

ji 


= lim sup — P 


n—)-co ^ 


1 


-(Tj - jmn) < yjjnin log( —) + 

V.7"in J 


i=i 


7 + 

Jn 


< limsup—= y(l + e). 

n^oo 


-5, 




By taking e J, 0 to the above and (|2.39l) we get (|2.29l) . The proof is completed. 


2.4 Proof of Theorem [3] 

Let {sj^r, lEj m} be independent random variables and sj^r ~ Beta(j, Z^)) 

that is, Sj^r has the Beta density gj,; < y < 1) for each y and r. Dehne 

m 

= n 1 < i < n. (2.40) 

r=l 

Let [x] denote the integer part of x and “A 0” indicate that “converges to zero in proba¬ 
bility”. 

We start with an auxiliary result before proving Theorem [3l 
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(2.41) 


Lemma 2.2 Let Gj{x) = P(Xn,j — ^ ^ [Oj l] <^'^d l < j <n. Then, 

Gi{x) > G2 {x) > ■ ■ ■ > Gn{x), X G [0,1]. 


Further, assume the conditions in Theorem 0 hold. If Ill.lM or lll.lll) is true, then 



(2.42) 


as n ^ oo, where and are as in Theorem\^ 


Proof. We first prove (j2.41l) . Let Xi and Yi be independent positive random variables for 


i = 1, 2. If P{Xi < x) > P{X 2 < x) and P{Yi < x) > P{Y 2 < x) for all x > 0, then one 
can easily show that P{XiYi < x) > P{X 2 Y 2 < x) for all x > 0. Combining this fact and 
(j2.40l) . it suffices to show that for each 1 < r < m and x > 0, P(sj,r < x) is non-increasing 
in j, where Sj^r is as in ()2.40p . 

Let Vj, j > 1 he i.i.d. random variables uniformly distributed over (0,1). For each 
n > 1, let Li:n < V2:n < ''' < ln:n denote the order statistics of Vi, • • • , 14,- It is well- 
known that Vj:n has a Beta(j, n — j + 1) distribution, see e.g., Balakrishnan and Cohen 
(1991). 

It is easy to see that Vi-n < V 2 :n-i-i < • • • < Vj.,n+j-i for any positive integers n and 
j, which implies P{Vj-n+j-i < x) is non-increasing in j for any positive integer n and 
X G (0,1). Since Beta(j, Z^) and have the same distribution, we have, for each r 

and X G (0,1), P(Beta(j, Z^) < x) is non-increasing in j. This concludes (I2.41h . 

Now we prove (I2.42p . Under condition lim^^oo I4(a^) = F{x) for x G [0,1] together 
with (|1.10l) or (|l.lll) . we first claim that 



(2.43) 


as n —)• oo for any x G (0,1). First, from (|2.30l) we have for any 6 G (0,1) 


-^^(1 “ ^) ^ io§^ ^ “(1 ~'t) for 5 < t < 1. 


(2.44) 


25 


Recall that Fn{x) = ( Y\^=i for 0 < x < 1. By assumption lim^^oo Pnix) = F{x) 

for X G [0,1], 



Since x < < 1 for any 1 < r < m, combining (I2.44p and the above equation we have 

that for any 0 < x < 1, 
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which implies (I2.43P since 1 — x > 0. 
Note that 


— All’ — f A 


j^r 


(2.45) 


j + Ir '' ' (j + lrY{j + /r + 1) 

Recall (j2.40p . From (jl.l9p and (j2.24p we can rewrite the moment generating function of 
log Ylj as 


M,(t) = H 


r(j + ^) r(j + ^r) 


L^iV r(j) r(j + t+ /,);• 


Hence we obtain 


77 m 

EdogFi) = = [M,d)^(logM,d))] = E ( 

r=l 


-r^(j) T^U + lr) 

*=o r(i + g 


that is, 


EilogYl^)=J2{m-Hj + lr)), 


(2.46) 


r=l 


where i/’(x) = r'(x)/r(x) is the Digamma function as mentioned in (I2.32p . Recall ?/(x) = 
X — 1 — log X for X > 0. Then 


log ^ log f^j,r = log ^ - 1) - ^ ^ ■ 


r=l 


From now on, for each x G (0,1), we take j = [nx]. We will show 


^ m 
1 


-E(—-1)^0 

7n 


and 


m 

7n /^i,r 


(2.47) 


(2.48) 


(2.49) 


To prove (I2.48p . it suffices to show that the variance of the left-hand side in (|2.48l) 
converges to zero. In fact, we have from p2.45p and p2.43l) 


Var( —V(^-l 
Wn ^ >i,r 


< 


1 0-2 

_L Ojx 

^2 ,,2 

'n r=l 

1 y. Ir 

7n ^ J(i + + 1) 

_ Ir _ 

7n ^ [nx]([nx] + Z,. + 1) 


0 ( 1 ) 

re7n f^^nx + lr 


E 
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as n ^ oo by the assumption 7n > 1- 

Since vi-r^) > 0, to show (|2.49p . it suffices to verify that t^E V™ , —S' 0 

as n —)• oo. To this end, we have from ()2.47p . p2.45p and (I2.46P that 



m 

( log - E log(y„^j) 

r=l 

m 

Y + ^r)) - (V’(j) - logj) 

r=l 



By formula 6.4.12 from Abramowitz and Stegun (1972) 


as f —>• + 00 . This and the fact j = \nx\ lead to 


^ m 

-eE^(—) 

7n ^ 


1 + 0 ( 1 ) 


In 


Y 


r=l •'J 


n'+^r- 1 


dt 


1 + o(l) l/i _ ^ \ 

In ^2 j' j+lr 

1 ±E 1 

7n 2 j(j + /r) 

m , 


7 n ^ nx + Zr 

r=l 


by (iTiSl) . 

From p2.47p - (I2.49P we have for any 0 < x < 1 


1 

In 



m 

^ ^ log+[nx],r^ ^ O' 

r=l 


Under pi.lOp or pi.lip , the limit F{x) is continuous and positive in (0,1). Therefore, the 
convergence lim„^ooEn(x) = F{x) is uniform for any interval [(5i,52] C (0,1). It follows 
that lim^^oo F„(il^) = F{x) for any x € (0,1). Further, notice that bn = 5n(l) and 
YT=i log+[nx],r = loggn(^) by (I2.45p . Then we have from ([OI) that 


J_ n,[nx\ _ 

'^n 


1 / ™ \ 1 o 

- (log +M - E log l-M-O + 7 ‘OS - ‘“S 

— ( log Yl[nx] - Y +[nx],r) + log - log F(x) 

^ r=l 
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converges in probability to zero for any x € (0,1). This completes the proof of (12.421) . ■ 

Proof of Theorem [3j Easily, part (c) is a corollary of (a) and (b). So we only need to 
prove (a) and (b). 

Notice, with the transform Z = iZe*®, that the density of Z is 2 ;^/*(kl) for 0 < |z| < 1 
is equivalent to that the density of (0,i?) is ^f*{r) for 9 G [0,27r) and 0 < r < 1. 

Let {sj^r, i < j ^ n,l < r < m} be independent random variables and sj^r have the 
Beta density < 2/ < 1) for each j and r. By Lemma [L3l for ease 

of notation we assume, without loss of generality, that 


Zj\ — 1 ^ J ^ 

r=l 


(2.50) 


If condition (jl.lOp holds, since Fn{0) = 0, Fn{l) = 1 and lim^^oo-^n(3:) = F{x) for 
X G (0,1), we assume without loss of generality that lim„^oo Fn{x) = F{x) for x G [0,1], 
F(0) = 0, E(l) = 1, and F{x) is continuous and strictly increasing over [0,1]. Hence, 
F*{y) := F~^[y) is a continuous and strictly increasing distribution function on [0,1] with 
F*(0) = 0 and F*(l) = 1. Define F*{y) = 1 for y > 1. Further, Unif{\z\ = 1} is identical to 
the product measure of Unif[0,2Tr) and hi. The distribution function of hi is the indicator 
function I[i^oo){y)- Let Yi, • • • , 1)1 be as in Lemma [LT] with ip = Wm’ defined in (|1.7p . 
According to Theorem [1] and Lemma l2.ll to prove (a) and (b), it suffices to verify that 


By Lemma ll.il 



Yly/-,. 

hn! 



F*{y) for all y > 0 if (11.1011 holds; 

/[I oo)(2/) for all y G (0,1) U (l,oo) if (11.111) holds. 



hn) 




(2.51) 


for any y G M. Since \Zj\‘^ G [0,1] is continuous for each j, it suffices to show 
1 ^ / 1 \Z-\'^ \ 

lim — P{ — log — < logy ) = the righ hand side of (12.511) . (2.52) 

n^oon^ ^7n K > 


From (I2.4ip and (|2.50p . 


( Y IZ'I^ \ 

P( — log — < log y) is non-increasing in j G {1, • • • , n}. 

V7n On ^ 

This property and equation (I2.42D play a central role in the following estimation. 


(2.53) 
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Proof of 112. 5 2\) as Hl.lll) holds. Fix y £ (0,1). For any 5 £ (0,1), we have 


1 

lim sup — P 

71^00 

[nS\ 




— log , 

In 0. 


< logy 


lim sup — r p (— log 

n-s.oo \jn 


i=i 


'f!l<iog!/)+ E p 

” i=M+i 


1 . I^7f 1 
— log < log y 

7n On 


< lim sup 


[n6] 


n^oo ^ 


+ lim sup p( — log 

n—>-cxD V 7n 




< 


log y) 


= 5, 




by the fact log A 0 as n —)■ oo for all x £ (0,1) from (12.421) and (II.lip . This implies 


^ ^ \^.\2 s 

lim — p( — log / < log y ) = 0 

n^oo Tl \ ’ 


i=l 


7n 


(2.54) 


for y £ (0,1). If y > 1, then logy > 0. By (12.5311 . 

1 > lim inf — y^ P{ — log —p— < log y 

rn-oo V7r,. On 

J=1 

^ [nr] 

> lim inf — y^ P 

n—>-oo n ^ ^ 


i=i 


nr 


1 , |z,-p , 

— log < log y 


> lim inf— -P [—log 


n^oo fi 


1 


\Z[n 


< 


In 


log y) 


= r 


as n —)• oo for all T £ (0,1). Letting r f 1, and combining with ()2.54p . we get (I2.52p . 

Proof of 112. 5 2\} as Hl.KX holds. We will differentiate two cases: y £ (0,1) and y > 1. 

Case 1: y £ (0,1). Set F*{y) = F~^{y) for y £ (0,1). Let 6 £ (0,1) be a number such 
that 0<y — 5<y + 5<l. Then 0 < F*{y — 6) < F*{y + 6) < 1, and F{F*{y — 6)) < y < 
F{F*{y + (5)). By letting x = F*{y + 6), we have 


1 


lim sup — y^ P 

n^oo 'O' . 

J=1 


1 

In 


■log 




lim sup [— y^ p{—\og < logy) + — y^ P 

n^oo In ^ \jn On 


< logy 


n—>oo Ln . 

J = 1 


j=[nx\ + l 


1 , |Zd2 

— log < log y 

‘^n bji 


< lim sup + lim sup P(— log < Jog y'j 

n—^oo ^ n—^oo '^7rx ' 


= x + hmsupP(^—log 

n^oo 

= F*{y + 6) 


\Z[n 


— log F(x) < — log iTT.'j 

V / 


(2.55) 
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by (I2.42P and the fact \og{F{x)/y) > 0. Similarly, setting x = F*{y — S), we get 


> 


> 


1 A /I |Z,f 

liminf — > Pi — log — < log 
n^oo n ^^ ^ h- 


J=1 

Inx] 


V7n 


S logy 


/I |Z,-M 

lim inf — > PI — log — < log y 

n^oo n \^n bn 


i=i 

\nx\ _ / 1 _ \Zu 


liminf 1 —-^P("—log ■ 

^In bn 


n—^oo Tl 

X 


< logy 


• lim inf P log -log F{x) < log 

n^oo V 7 „ bn F{x)J 


X ■ 


-'n 

1 \Zu 


1 — lim sup P ( — log ■ 

n^oo '^7n 


P(a 

— logP( 3 :) > log 


F{x) 


= F*{y-6) (2.56) 

by (I2.42P and the assertion log(y/P(x)) > 0. Finally, by letting 5 | 0 in (12.551) and (12.561) . 
we show (12.521) holds under the condition (II.IUI) and y G (0,1). 

Case 2 : y > 1. Observe 


1 |'7|2 1 

-^pf—log < logy) > -^p(—log 

Vyn bn 


\Z, 


< log yi 


for all yi G (0,1). By the proved conclusion, it is seen that 


1 ^ / 1 \Z\ 

liminf — P( — log —-— < logy ) > P*(yi). 

n^oo \jn bn > 

J=1 

Then (j2.52l) follows by taking yi t 1 a'lib by the fact P*(l) = 1. The proof is complete. ■ 

Now we present the proofs of the corollaries. 

Proof of Corollary [TJ Since Uj = Ij + n, take 7 ^ = 2 to have 


^n(x)=(n 


rijX 


Ui — nil — x) 
1=1 ■' 


1/2 


0 < X < 1. 


By assumption, lim^^oo ^ [Oj 1] for 1 < j < m, we have 

/Th X \V2 

lim P„(x) = ( TT ----=: P(x), 0 < x < 1. (2.57) 

n^oo — 0,(1 —X)/ 

j=l J 

(1). If oi = • • • = am = 1, then P(x) = 1 for 0 < x < 1. That is, condition (|l.lll) holds. 
The conclusion follows from (c) of Theorem [3l 
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(2) & (3). If Uj < 1 for some 1 < j < m, then F{x) in (I2.57P is continuous and strictly 
increasing in (0,1]. Also, lima;j_o = 0. The statement (3) then follows from (a) of 
Theorem [3l In particular, when ai = ■ ■ ■ = am = a G [0,1), then F{x) = ( 
for 0 < X < 1. It is trivial to check that 


F-\x) 


(1 — a)x^/'^ 

1 — ax‘^1'^ 


and 


±F-Ur\ = 2 ( 1 -a) 

dx m [1 — ax^l'^Y 


We get (2). ■ 

Proof of Corollary [2j Take 7 ^ = m,„. We will show that lim^^oo = F{x) with 


F(x) = xexp^ —y log{l — q{t){l — x))dt^, 0 < x < 1, (2.58) 

and lim 2 ; 4 ^o-^( 21 ) = 0. Obviously, T(l) = 1. This says that condition (jl.lOp is satisfied. We 
first prove (I2.58p . Since 


logF„(x) 


1 >'‘'71 

—5:[iogx-iog(i-^(i-x)) 

j=l J 

logx-Vlog (1 - (1 - x)), 

rrin^^ V Uj J 

1=\ 


we have that 


logFn(x) - logT(x)| = 


< 


1 / \ 

— (1 - —( 1 -a;)) - / log(l - g(t)(l - x))dt 

fTlfi \ Tlj / Jq 

-| rUn 

7; E [l»8(l - 7(1 - *» - ‘“8 (1 - ‘l(i) (1 - *)). 

' 1=1 


^ TTLn . 

— ^ log (1 - q(^) (1 - ^)) - / log(l - q{t){l - x))dt 

11 In ■ -I 8 \iiLyi' ' JQ 


+ 

777 ,„ 

i=i 

: /i(x) + hix) 


for 0 < X < 1. Since for each fixed x G (0,1], ^ log(l — t(l — x)) = for 0 < t < 1, 

we have 1 — x < log(l — t(l — x)) < By using the mean-value theorem and (I1.12p . 
we have 


h{x) < 


1 — X 1 


X 


iitn 


i=i 


\mr 


0 


as n —)• CO. For fixed x G (0,1], since log(l — t(l — x)) is a bounded and continuous 
function in t G [0,1] and thus log(l — q{t){l — x)) is also a bounded and continuous function 
in t G [0,1], we have Fix) —)• 0 as n —)• 00 by the definition of the Riemann integral. 
Therefore, we have proved that limn^.oo Fnix) = F{x) for x G (0,1]. 
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Choose any to € (0,1) such that (?(to) < 1- From continuity, there exists e > 0 satisfying 
0<to — €<to + e<l and sup|^_^|j|<g q{t) < 1 — e. Easily, — g(tKi-x) < 1 for all X G (0,1). 
It follows that 


logF(x) = logx— / log(l — g(t)(l — 




X 


-dt 


-q{t){l-x) 

/■io+e 

Jto-e - x) 


X 


dt 


< log 


2ex 




—oo 


as X I 0, yielding that lim^-j^o -F(x) = 0. Thus, (ll.lOp is verified. 
Now, 


/(x) = F'(x) = ^F{x)-F{x)[\ 
X Jo ^ 


qit) 


0 

p, s, 1 -^(i) 

I 


dt 


Further, 


d 


rw = = 


0 l-9(i)(l-a^) 
1 


dt. 


f{F-Hx)y 


0 < X < 1. 


By (a) of Theorem[3l d-j where fj, has density ^f*{r) for {6,r) G [0,27r) x [0,1). 

Proof of Corollary [3j Fix 0 < x < 1. Then, 


1 m-n ^ ^ 

logF„(x) =-y^log -(iH-(x-1 

7n ^ ^X\ no 


7n— LX 


(2.59) 


Write 


n 


1 + —(x-1) = 1 + 


1 — X 


1 - 


n 


n. 


By the assumption lim„^ooiiiaxi<j<m„ |^ — 1| = 0, we know the logarithm on the right 
side of (|2.59p is equal to ^^(l — ^)(1 + o(l)) uniformly for all 1 < j < mn- Use the 
identity that rij = lj + n to see 

1 TJ'i'n 1 

logF„(x) = 3^(1 + o(l)) ^ ^ 


InX 


rij 


1=1 


, , ,,x-l 1 1 

= (1+ o(i))—• — • 0 < X < 1, 


X 7n n 


(2.60) 


1=1 


35 


















as n —)■ oo by the assumption lim^^oo |^ — 1| =0 again. We now show (a), (b) 

and (c) according to the assumption that lim^^oo ^ ET=iij = P G [ 0 ,oo]. 

(a) . Assume /3 = 0. Choose jn = 2. Then, T’(x) = 1 for 0 < x < 1. The conclusion follows 
from the second part of (c) from Theorem [3l 

(b) . Assume /3 £ (0,oo). Choose 7 „ = 2. Then F(x) = exp(^^^) for 0 < x < 1 and 
liniajj^o T(x) = 0. Trivially, F~^(x) = (l — ^ logx) ^ for 0 < x < 1 and 


d 


F-Hx) = 


2/3 


dx^ x(/3 — 21ogx)2 ’ 

So the density of /i according to (a) of Theorem [3] is 

/3 


0 < X < 1. 


1(0,r) := 


7 rr (/3 — 2 log r)^ 


for {6,r) G [0,27r) x [0,1). By the polar transformation z = re*®, it is easy to see that the 
density of fi is given by 

H^) ■= I lo/a ^oi - Ti^’ 1^1 < 1- 

7 r| 2 :p(p — 21 og \z\)^ 

(c). Assume /3 = oo. Take 7 n = 7 h- Then, by (|2.60p . F{x) = exp{^^} for 0 < x < 1 
and lima; 4 ,o F{x) = 0. This is the same as the case that /3 = 2 in (b). So the density of p is 
„ I —nw for 0 < jz] < 1. ■ 

Proof of Corollary [4l By assumption, m = — )• 00 and maxi<j<m„ ^ = 0 as n —)> 00 . 

By using the equality Uj = Ij + n and taking 7 ^ = mn, we see that 


lim Fn{x) = X := F{x) 
n^oo 

for any x G [0,1]. Then F~^{x) = x for x G [0,1] and /*(x) = 1 for 0 < x < 1. Choosing 
hn{x) = ajid 5 ^ = YYP=i ■'^0 know that fj-n and Z = iZe*® has 

the density 2 ^]^ for 0 < jzj < 1. Equivalently, (0,i?) has law 17m/([0,27r) x [0,1]). As 
mentioned before, Z = has the uniform distribution on {z G C; \z\ < 1} if and only if 
(0, \Z\'^) follows the uniform distribution on [0, 2it) X [0,1). Therefore, the desired conclu¬ 
sion follows. ■ 
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